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ABSTRACT 


A  variant  of  the  simplex  method  Is  given  for  solving  linear  pro¬ 
grams  with  M  +  L  equations,  L  of  which  have  the  property  that  each 
variable  has  at  most  one  nonzero  coefficient.  Special  cases  Include 
transportation  problems,  programs  with  upper  bounded  variables,  as¬ 
signment  and  weighted  distribution  problems.  The  main  feature  of  the 
algorithm  is  that  a  working  basis  of  order  2M  +  1  is  utilized  for 
pivoting,  pricing,  and  inversion  which  for  large  L  can  be  of  signifi¬ 
cantly  lower  order  than  that  of  the  original  system.  A  method  utilizing 
a  working  basis  of  order  M  in  which  the  updating  of  the  inverse  re¬ 
quires  elementary  rcw  as  well  as  column  transformations  will  be  dis¬ 
cussed  in  a  sequel  ^  . 


GENERALIZED  UPPER  BOUNDED  TECHNIQUES  FOR  LINEAR  PROGRAMMING  -  I 


I.  Introduction 

We  are  interested  in  linear  programs  with  M  +  L  equations  such 
that  each  variable  has  at  most  one  non- zero  coefficient  in  the  last  L 
equations  [Figure  1]  .  Transportation  problems,  linear  programs  with 
upper  bounded  variables,  assignment  and  distribution  problems  all  have 
this  structure.  To  simplify  the  discussion,  we  assume  all  the  coeffi- 
cients  in  the  last  L  equations  to  be  non-negative.  The  constant  terms 
are  assumed  to  be  positive.  In  section  6  the  necessary  modifications 
to  handle  negative  coefficients  in  the  last  L  equations  will  be 
outlined. 

By  dividing  each  of  the  last  L  equations  by  the  corresponding 
right-hand  side  element  and  by  scaling  the  variables  in  the  problem 
we  can  assume  without  loss  of  generality  that  each  equation  in  the 
subset  is  a  partial  sum  of  the  variables  equal  to  unity.  See  figure  2. 

In  the  next  section,  some  terminology  is  introduced  and  a  few 
theorems  given  which  are  the  basis  for  the  proposed  algorithm.  The 
method  to  be  used  is  explained  in  the  third  section  and  in  the  following 
section  the  algorithm  is  set  forth  in  detail.  In  the  fifth  section 
the  algorithm  is  illustrated  by  an  example.  In  the  final  section  the 
slight  modifications  needed  to  handle  negative  coefficients  in  the  last 
L  equations  are  given. 

II.  Some  Theorems  and  Definitions 

th 

The  i  set  of  variables  or  columns.  S^  ,  will  refer  (depending 
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subject  to 


on  context)  to  those  variables  or  columns  corresponding  to  the  columns 

th 

of  coefficients  in  figure  2,  with  1  as  their  M+i  component.  SQ  , 
th 

the  0  net .  is  the  set  correxponding  to  columns  with  zeros  for  the 

st  th 

m+1  through  M+L  coefficients. 

We  assume  that  the  system  in  figure  2  is  of  full  rank  and  denote 
1  M+L 

by  [B  , ...,B  ]  a  basis  for  the  system.  The  underscoring  is  to 

differentiate  coefficient  vectors  with  all  M+L  of  their  components, 
from  reduced  vectors  to  be  introduced  shortly  which  are  2M+1 -vectors. 

There  will  be  no  underscoring  for  individual  coefficients  since 
the  two  different  types  of  vectors  differ  only  in  the  number  of  components. 
THEOREM  l)  At  least  one  variable  from  each  set  3^  (i  ■  1,...,L  ; 
excluding  SQ)  is  basic. 

PROOF:  £31  nee  we  have  assumed  that  our  system  is  of  full  rank  any 

M+L-vector  can  be  represented  as  a  linear  combination  of  vectors  in  the 

basis.  Consider  the  vector  of  ones.  13iere  are  X,  such  that 

J  j 

I  -  [l,...,l]T  .  In  order  for  2  *  1  >  B^  must  be 

J 

nonzero  for  some  and  this  can  only  happen  if  B  €  . 

THEOREM  2)  The  number  of  sets  containing  two  or  more  basic  variables 
is  at  most  M  . 


PROOF:  By  the  assumption  of  full  rank,  each  basis  has  exactly  M+L 
vectors.  By  Theorem  1,  L  of  them  are  in  separate  sets,  this  leaves 
at  most  M  to  make  up  sets  of  more  than  one  basic  variable. 

According  to  Theorem  2  we  can  partition  all  the  sets  excluding 
SQ  into  those  containing  exactly  one  basic  variable  and  those  containing 
more  than  one  basic  variable.  At  any  given  point  in  the  algorithm  we 


r  * 


will  call  M+2  of  the  sets  essential.^  Included  In  these  sets  are  all 
the  sets  with  more  than  one  basic  variable  plus  the  set  .  The 
remaining  sets  are  called  inessential. 

REMARK:  By  Theorem  2,  at  least  one  of  the  essential  sets  besides  ^ 
has  exactly  one  basic  variable  in  it. 

Columns  of  coefficients  and  variables  are  called  essential  (e£) 
or  inessential  (^£)  in  a  similar  manner. 

We  define  the  reduced  system  to  be  the  system  obtained  from  the 
system  in  figure  2  by  deleting  the  equations  corresponding  to  the  rows 
of  "1"  coefficients  for  the  Inessential  sets.  The  set  of  essential 
basic  columns  restricted  to  the  reduced  system  is  called  the  working 
basis.  Inessential  basic  variables  are  said  to  be  at  their  upper  bound 
since  their  value  must  be  one  for  any  basic  feasible  solution.  Columns 
restricted  to  the  reduced  system  are  denoted  by  to  distinguish 
them  from  the  which  is  the  whole  column  of  MfL  components. 

Since  ve  had  M+L  equations  initially  and  removed  L+l-  (M+2)  equations, 
A"  has  2M+1  components. 

THEOREM  3)  The  working  basis  is  a  basis  for  the  reduced  system. 

PROOF:  The  number  of  equations  in  the  reduced  system  equals  the  number 

of  variables  in  the  working  basis  since  for  each  inessential  set  there 

is  exactly  one  equation  removed  to  form  the  reduced  system  and  exactly 

one  inessential  basic  variable  which  is  not  Included  in  the  working 

basis.  All  that  remains  is  to  show  that  the  columns  of  the  working 

basis  are  linearly  independent.  Suppose  not,  then  there  exists, 

1  We  make  the  assumption  that  L  >  M+2  and  as  a  matter  of  fact  the 
algorithm  will  only  be  practice],  if  L  »  M+2  . 
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I 


^  *11  zero,  such  that 


z 


■  0  where  (B^)  J  ■  1,...,2M*1 


are  the  columns  of  the  working  basis.  But  each  B^  in  the  reduced 
system  differs  from  BJ  ,  the  whole  column,  only  by  zero  components  so 
that 


I 


¥ 


J 


contradicting  the  linear  independence  of  the  original  basis. 

The  three  theorems  proved  In  this  section  allow  us  to  work  with 
a  reduced  system  of  2M+1  equations.  If  L  »  Mf2  this  will  allow 
us  to  do  our  simplex  operations  In  a  system  of  a  much  smaller  order  than 
the  original. 


III.  The  Method 


J 


Suppose  we  have  a  basis  B  ■  (B  :  J  ■  l,...,MfL}  for  the  whole 
system  where  we  assume  without  loss  of  generality  that  B1  ■  Ap  ,  the 
column  corresponding  to  the  variable,  x^  ,  to  be  maximized.  Let 

ki 

B  -  (B  J  -  1,...,2M+1) 

k  k  k 

be  the  working  basis  where  B  ”  c  B  and  B  ^  is  obtained  from  B  ^ 

kl  l  o 

by  deleting  the  appropriate  components.  Again  we  assume  B  -  o  ■  A  . 
Oir  first  project  is  to  find  the  value  of  the  basic  variables  for  the 
basic  solution  determined  by  B.  That  is  we  seek  x^  k  ■  1, . ..,M+L 
such  that  L  -t*  Let  U  denote  the  set  of  variables  at  upper 
bound,  l.e.  inessential  basic  variables.  Then  we  have 
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+  X*/-  5V + 

B*Vu  B^eU  B^/U  B^eU 


since  we  know  that  x^  ■  1  for  B^  c  U  ,  all  that  remains  is  to  find 

the  value  for  variables  in  the  working  basis  i.e.,  x.  such  that 
J.  Jk 

B  K  /  U  .  Consider 


The  columns  on  the  left-hand  side  of  the  above  equality  are  nonzero  at 

most  in  the  components  corresponding  to  the  reduced  system.  Similarly 

on  the  right  side,  for  S.  which  are  inessential  b^  •  1  and  there 

1  k  k  “Mfi 

exists  exactly  one  column  B  0  c  U  with  ■  1  ,  Thus  to  find  W  , 

the  values  of  the  variables  corresponding  to  the  working  basis,  it 
suffices  to  consider  only  the  reduced  system,  i.e., 


w  -  [ \  J  -  B  |b  -  Y  bM 

L  *1  *2M+1J  L  V  J 

r€u 


The  next  item  of  interest  is  to  find  the  prices  corresponding 
to  the  current  basis.  Let  denote  the  prices  on  the  first 

m  equations  and  ^,...,1^  the  prices  on  the  last  L  .  Since  we 
assumed  that  B1  ■  fP  the  first  row  of  B  ^  ,*  denoted  (B1)^  ,  will 
be  a  set  of  prices  for  the  reduced  system  in  the  sense  that  (B’^)^B  ■ 
[1,0,0, ...,0]  a  2M+1  conr±x>nent  vector.  We  now  extend  these  prices 
to  a  set  of  prices  for  the  whole  system.  Let  7Tj  •  (B  j  •  1,...,M 
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Now  let  us  see  what  the  i  ■  1,...,L  have  to  be.  If  3^  is 

essential  and  e  then 

M±L 

o .  x  (b*1); 

1-1 

2Mfl 

■  l  (B_1)J  Bi 

4-1 

- 1  <.->>; .  i 

i-i 

■  L  "A  +  "I  * 

Hence  ^ 

(2)  ^  f0r  -  €  ^ 

i-1 

The  jjl  ,  so  obtained,  are  simply  the  appropriate  components  of 
(B  X)1  .  However,  the  same  formula  (2)  can  be  used  to  extend  our 
prices  to  the  full  system  Including  inessential  sets.  This  is  clear 
because  the  prices  [7r,(i]  are  determined  by  the  requirement  that 

]_  i  i 

B  *  1  ,  and  [7T,m.]B  »  0  J  -  2,...,M+L  .  For  2°  e  £  these 

equations  hold  because  n  was  defined  in  terms  of  B  ^  .  For  £ 
the  requirements  are  Just  (2)  again.  Since  every  inessential  basic 
variable  appears  in  exactly  one  equation  of  the  form  (?)  [7r,pJ 
is  a  set  of  prices  for  the  original  system. 

Now  we  are  able  to  "price  out"  the  columns  of  the  original  system 


-8- 


to  see  which  one  should  come  into  the  basis,  e.g.,  using  the  usual 

s 

simplex  criterion  we  might  choose  the  Incoming  column  A  by 


A 

s 


(7T,h)AS 


min  (tt.u)a^ 

J 


min  A. 

J  J 


If  As  £  0  we  have  an  optimal  basic  feasible  solution,  if  not  we  bring 
s 

A  into  the  basis.  We  note  that  -he  ii^  need  not  be  stored  in  the 
computer  but  may  be  generated  as  needed  by  (2).  Let  us  suppose  that 

*  V 

If  is  essential  we  simply  pivot  in  the  reduced  system  using 
the  usual  criterion  for  picking  the  row,  i.e.,  the  row,  r  ,  is  given  by 


(b’V^x) 

where  of  course  (B  A  )  >  0  ,  and  we  only  consider  rows  in  the  reduced 

system,  except  r  *  1  since  the  column  corresponding  to  the  functional 
always  remains  in  the  basis.  We  then  proceed  with  another  iteration. 

If  is  inessential  two  things  can  happen : 

(1)  If  AS  drives  out  of  the  basis,  where  B^  is  the 

g 

unique  basic  variable  in  then  A  simply  goes  to  upper 

t 

bound  instead  of  B  and  we  need  not  pivot  at  all.  Just 
modify  the  sum  b  -  I  B^  in  (l)  and  (4). 

(2)  If,  on  the  other  hand,  A8  drives  out  some  other  basic 

(and  necessarily  essential)  variable,  we  redefine  the  essential 


[B'V  -  I  BJ)] 

B'U  r 

(B'V) 

'r 


©  ■  min 


[B'V, I  BJ)] 
BJcU  1 


(BA8), 
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sets  so  as  to  make  S  essential  and  then  perform  a 

0 

pivot  using  (4). 

Nov  let  us  go  into  this  in  more  detail.  In  particular,  to  see 

whether  (l)  or  (2)  occurs, ve  let  x  ■  0  be  the  Increase  in  the  variable 

x  .  Simultaneously  by  the  sum  constraint  x  ,  the  coefficient  of 
s  t 

,  is  decreasing  from  x  *  1  to  x  ■  l  -  Q  ,  To  satisfy  the 

w  w 

equations  in  the  reduced  system  ve  must  adjust  the  variables  of  the 
working  basis.  Equation  (l)  now  becomes 


«(e) 


-  b'1  b  -  V  -  8(A®  -  Bk) 


„k 
B  cU 


(5) 


-  B+  -  ^  fl*]  -  OB*1  AS  -  B*] 
^eU 


«  W(o)  -  OB 


-1 


A  -  B 


which  is  required  to  be  non-negative  except  possibly  in  the  first  com¬ 
ponent  which  corresponds  to  the  functional.  If  ve  let  the  i-th  component 
of  W(o)  be  Wi  ,  then  the  maximum  0  for  feasibility  Is 


(6) 


0  ■  min 

tfi 

[B’V-BbL  >  o 


[B'V-B1)]. 


[B'1(AS-Bt)]. 


where  IB  ^(A®  -  B1)]  >  0  • 
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Hence  0  a  min [9,  1]  where  W(Q)*  (x.  )  are  the  values 

K1  K2M+1 

for  the  variables  in  the  working  basis  from  the  previous  iteration. 

Of  course,  the  minimum  only  makes  sense  over  the  rows  in  the  reduced 

system.  If  9  »  1  then  x  goes  to  upper  bound,  and  the  variable 

s 

t 

corresponding  to  B  goes  to  zero.  If  we  indicate  by  a  super  bar  the 
new  values  cf  the  various  variables  we  obtain: 


w  r,  w  -  b'V2  -  B^ 


If  on  the  other  hand  9  <  1  we  first  make  S  essential  (which 

o 

4- 

automatically  introduces  B*  into  the  working  basis)  and  then  introduce 
s 

A  into  the  new  working  basis.  To  do  this  we  must  make  one  of  the 

currently  essential  sets  inessential.  This  essential  set  cannot  be 

SQ  and  clearly  must  contain  exactly  one  basic  variable.  Put  by  the  remark 

maae  after  theorem  2  there  must  exist  such  a  set.  Denote  by  S  one 

k  P 

I* 

such  set  and  let  B  be  the  basic  variable  in  S  . 

P 

— 1 

Uow  we  must  obtain  an  inverse  E  for  the  new  working  basis, 

-  -  -1 

B  .  B  can  be  defined  by  (B  )  B  J  ■  5  where  6.  ,  ■ 

^  1J  J 

0,1  J  and  5  a  1,  i  ■  J  ,  For  the  old  inverse  B  we  have 

1/ 

-1  K 1  — 

(B  )  B  j  •*  6.  ,  .  But  B  differs  from  B  only  by  the  fact  that 

1  _  1J  k 

t 

column  B  replaces  B  and  the  row  M  ♦  a  replaces  M  +  p  in 

t  k 

the  reduced  system.  But  except  for  B  and  B  r  all  the  columns 
involved  have  0's  in  both  the  M  +  p  and  M  +  0  th  components. 
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Let  denote  the  matrix  obtained  by  the  permutation  of  the  columns 

of  B~^  corresponding  to  the  changes  in  the  rows  of  the  working  basis 

th  th 

resulting  from  replacing  the  M  ♦  p  equation  by  the  M  ♦  o  in  the 

-1  -ki 

reduced  system.  We  have  then  (B  )'  B  d  -  6.  for  all  k.  ^  t  .  To  get 

11  1 J  J 

—  •! 

B  all  ve  have  to  do  is  pivot  on  B  )  'B  on  the  element  corresponding 

to  the  M  ♦  o  rcw.  This  pivot  is  particularly  easy  since  the  pivot 

element  will  always  be  one,  which  can  easily  be  verified.  U  the  set 

k 

of  variables  at  "upper  bound"  is  now  changed  by  including  B  r  and 
t  - 

deleting  B  and  V  ,  of  course,  remains  unchanged,  except  possibly 

_-l 

for  the  order  of  the  components.  Using  the  new  B  we  could  obtain 
new  prices  which,  however, mist  be  the  same  as  the  old  ones  since  the 
basis  for  the  original  system  is  the  same.  Again  A6  will  price  out 
optimally  but  new  is  essential  so  we  Just  do  an  ordinary  pivot  on 

rew  r  determined  by  (6),  and  we  go  on  to  the  next  iteration.  This 
essentially  describes  the  method.  In  the  next  section  the  algorithm 
is  described  in  detail. 


IV.  Description  of  Algorithm: 

Referring  to  figure  the  algorithm  takes  place  in  the  following 
steps : 

(1)  We  assume  we  enter  the  algorithm  with  the  inverse  B-1  of  the 

working  basis  B  ,  the  value  of  the  variables  in  the  working 

basis  W  =  B  ^(b  -  X(B^)  and  the  set  of  variables  U  at  upper 

BJeU 

bound.  To  get  this  initial  solution  the  usual  phase  I  procedure 
can  be  used  in  the  obvious  way. 

(2)  Let  7^  =  (B**)*  for  i  «  1,...,M  where  we  assume  for  simplicity 
that  the  column  corresponding  to  the  objective  variable  is  the 
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first  column  in  the  working  basis.  For  each  set  S  (v  i  0)  let 

\a  a  a  V 

M  3  J 

Aj  ■  121  7 A"  +  ^  for  A0eSv  .  Let  A§  ■  min  Aj  and  suppose 
Afl  €  6^  .  If  Ag  £  0  we  go  to  step  (3),  if  A#  <  0  and  SQ 
is  essential  we  go  to  step  (4).  Finally  if  A  <  0  and  S 

0  Q 

inessential  ve  go  to  step  (6). 

(3)  Terminate  ve  have  an  optimal  solution. 

( h )  Ve  find  A6  in  terms  of  the  working  basis,  that  is  B*^  A8  and 
find  the  pivot  row,  r  ,  by  equation  (4). 

(3)  Ve  then  pivot  on  row  r  column  s  in  the  reduced  system  giving 

_  _-l 

us  a  new  working  basis,  B  ,  and  inverse,  B  .  After  updating 
by  (1),  ve  return  to  step  (1)  for  the  next  lnteratlon. 

“1  a  ♦  ^ 

(6)  Ve  find  B  (A  -  B  )  and  determine  the  value  of  9  and  row  r 

by  equation  (6).  If  0  >  1  then  x  *  0  ■  1  and  ve  go  to 

®  0 

A  A 

step  (7).  If  9  <  1  then  x  ■  9  ■  9  and  we  go  to  step  (8). 

0 

-8  t 

(7)  A  goes  to  upper  bound  in  set  and  B  the  unique  variable 

which  was  at  upper  bound  goes  to  zero.  The  working  basis  and  the 
reduced  system  remain  the  same  but  V  and  the  prices  are  updated 
by  equations  (7). 


7riAiV  vhere  V  a  ^ai^c  column  in  .  Let 


(b)  Suppose  S  (p  /  0)  is  an  essential  set  with  exactly  one  basic 

Pk  A  ^ 

variable  B  r  .  Let  bJ«b|ij<M  +  o,M  +  p  and  ^  ■  0 

^  t  t 

and  B„  *  1  .  That  is  ve  move  the  "one"  coefficient  in  B 
M+p  k 

to  the  same  row  as  the  one  coefficient  in  B  r  .  Then  ve  pivot 
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in  the  reduced  system  on  B 


in  the  M  ♦  p  row  regarding 


B  as  belonging  to  .  We  then  permute  the  columns  of  the 

resulting  inverse  so  that  we  can  replace  the  M  ♦  by 

th  *r 

M  +  o  equation  in  the  reduced  system  obtaining  B  .  B 

t 

replaces  B  in  U  and  the  elements  of  W  are  permuted  to 

_-l 

reflect  the  column  reordering  of  B 


Q 

(9)  New  we  need  to  find  A  in  terms  of  the  working  basis  obtained 

--1  8  --1 

in  step  (8),  that  is  we  want  B  A  .  If  B  is  obtained  from 

B_1  by  T  a  pivot  transformation  matrix  and  P  a  column  per- 

-_1  -1 

mutation  matrix  then  B  -  PTB  .  But  we  have  already  computed 
B  ^(4S  -  Bt)  in  step  (6)  and  PTB  Si*  =  U  a  unit  vector  with  a 
one  in  the  row  of  the  reduced  system  corresponding  to  the  M  ♦ 
equation  of  the  original  system,  by  the  way  P  and  T  were 
chosen,  hence: 


BA8-  B  (A®  -  B1)  ♦  U 
*  PTB’1(A2  -  At)  ♦  U 

-1 

—  6  /v 

we  then  pivot  on  B  A  as  in  step  (5;  using  the  row  r  deter¬ 
mined  in  step  (6)  . 
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ENTER 


Figure  J 
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V.  Example :  Consider  the  following  example  for  M  *  2  with  the  co¬ 


efficient  array  given  in  figure  4.  We  are  trying  to  maximize  x. 


A1  A2  A}  A4  A5 


A7  A^ 


The  initial  reduced  system  contains  variables  xQ  ,  x^  ,  x^  ,  x^ 

5  7 

and  the  first  five  equations.  U  =  {A  ,  A  } .  The  inverse  of 

0  12  5  4 

B  «  A  A  A  At  A  is 


4  4-4  -5/2 

-i  i  3/4  i 

4-4  4-4 

0  0  0  1 

0  0  0  0 


Fran  the  first  row  of 


ve  obtain 


n 


-5/2 

3/4 

-3/4 

0 

1 

[1/2,  1/2]  and  by 
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pricing  out  basic  columns  we  obtain  p.  «  (-1/2,  -5/2,  -5/2,  m)/2,  2]  , 
hence  [II,  \i]  =  [1/2,  1/2,  -1/2,  -5/2,  -5/2,  -9/2,  2]  is  the  pricing 
vector.  After  pricing  out  we  find  *  -3  is  the  minimum.  Since 
e  which  is  inessential  we  find 


B"1^6  -  A'’ ) 


-3 

1/2 

-1/2 

0 

0 


and  we  use  (6)  to  compute  ©  =  =  1  since  there  is  only  one  posi- 

tive  coefficient  in  B~  (A  -  A?)  .  Since  ©  £  1  ,  goes  to  upper 
bound  and  x<-  goes  to  zero.  The  new  values  for  the  variables  in  the 
working  base  are  given  by 

W  =  (x0,  Xx,  x2,  Xy  x/  =  h  -  B~X(A°  -  A5) 

=  [3,  1/2,  1/2,  1,  1]T  -  [-3,  1/2,  -1/2,  0,  0]T 
=  [6,  0,  1,  1,  1]T 

and 

x  =  (6,  0,  1,  1,  1,  0,  1,  10]T 

[ff,  £]  =  [1/2,  1/2,  -1/2,  -5/2,  -5/2,  -3/2,  2] 


since 


-  A6  *  -9/2  ♦  3  =  -3/2  . 
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Q 

New  we  start  over  with  the  new  prices.  This  time  A  €  prices  out 
optimally  with  Ag  *  *1  .  S,  is  also  inessential  so  we  compute 

B‘1(A8  -  A7)  -  [-1,  5,  -5,  o,  0]T  . 

Since  (B*^(A8  -  A7)^  is  >  0  only  for  1  ■  2  we  have  ©  ■  j  <  1  . 

Which  tells  us  we  have  to  make  essential.  To  do  this,  we  oust 

make  one  of  the  previously  essential  sets  inessential.  Sg  and 

are  both  available  for  this  purpose  since  they  both  contain  exactly 

one  basic  variable.  Let  us  move  Sg  out.  The  new  reduced  system 

will  include  the  first,  second,  third,  fifth,  and  seventh  equations. 

-‘1  -_1-j 

We  new  seek  the  uew  inverse  B  so  that  B  Au  is  a  unit  vector 

for  J  =  0,  1,  2,  4,  7  .  To  do  this  we  consider  A'^  *=  [-1,  -3,  0,  1, 

where  we  have  changed  the  one  from  the  seventh  component  to  the  foirth. 

7 

now  represent  A"  in  terms  of  the  current  basis: 

B'V7  =  [-9/2,  -1/4,  +l/4,@,  0]T 

and  pivot  on  the  fourth  component.  The  matrix  which  accomplishes  the 
pivot  is 


T  « 


1  9/2 

1  1/4 

1  -1/4 
1 

0  1 


and 


,  0,  0]T 

We 
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B 


1/2 

1/2 

-1/2 

2 

-5/2 

-1/4 

1/4 

3/4 

1/2 

3/4 

1/4 

-1/4 

1/4 

-1/2 

-3/4 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

In  order  to  use  the  above  as  the  inverse  working  basis  we  obtain 

from  A^  by  taking  A*^  =  [AjJ,  A^,  Ay  Ay  A^]T  however  if  we  inter- 

-1  — 

change  the  fourth  and  fifth  row's  and  columns  of  B'  ,  obtaining  B 


then  A 


J 


[a: 


A^,  A^,  A^,  A^]  i.e.,  in  ascending  order  of  the 

1  “1  A  “1  1  “1 


com¬ 


ponents  in  ,  and  the  product  B  A^  *  U.  ,  B  A^  *  U0  ,  B  f?  «  U,  , 
-1-1  x  2  3 

-  4  -  x  7 

B  A  •-*  and  B  A  =  ,  where  is  the  unit  vector  with  1  in 

its  1th  component  and  zeroes  elsewhere. 


1/2 

1/2 

-1/2 

-5/2 

2 

-1/4 

1/4 

3/4 

3/4 

1/2 

1/4 

-1/4 

1/4 

-3/4 

-1/2 

0 

0 

0 

1 

0 

0 

0 

0 

0 

1 

-1 


—  A  8  -1  8  7 

New  that.  S^-  is  essential  we  can  pi\ut  on  B  A°  *  PTB*  (A  -  A')  + 
where  P  is  the  permutation  matrix  that  performs  the  interchange  of 


the  fourth  and  fifth  columns  and  rows. 


B  AC  .  PT[-1,  5,  -5,  0,  0]T  +  Uj 
-  -5,  0,  1JT 
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The  values  of  X  are  the  sa me,  but  the  components  which  are  in  the 
working  basis  change;  i.e.. 

Previously  W  *  [xQ,  x^  Xg,  Xy  x^]  -  [6,  0,  1,  1,  1] 

Now  W  -  [xQ,  xx,  Xg,  x4,  x?]  -  [6,  0,  1,  1,  1]  . 

__-l  Q 

If  we  now  pivot  on  the  second  component  of  B  A°  using 


1 


+1/5 

1/5 

1  1 


o  l 


l 


and  interchange  the  second  and  fifth  r^s  by  using  P1  to  get  the  unit 

_-l  < 

vectors  B  Ad  ,  J  =  0,  2  ,  4,  7,  8  in  numerical  order  on  J  ;  we  obtain 


9/20 

11/20 

-7/20 

-47/20 

42/20 

_-l 

1/20 

-1/20 

-3/20 

-3/20 

18/20 

B  * 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

-1/20 

1/20 

3/20 

3/20 

2/20 

n  =  [9/20,  11/20] 

and 

4  *  [-7/20,  -47/20,  -47/20,  -31/20,  42/20] 
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The  full  pricing  vector  is  then 


[n,  m3  =  [9/20,  11/20,  -7/20,  -47/20,  -4 

Upon  pricing  out  we  find  L.  ^  0  for  all  j 

J 

optimal  solution,  the  value  of  which  we  obtain  by 


/20,  -31/20,  42/20]  . 


and  hence  have  an 


=  [-*0’  X2’  xl»*  X?*  X8^ 

= 

=  [6,  1,  1,  1,  0] 


and  x,  =  x,  =  1 
o  3 


are  at  upper  bound. 


So  the  full  optimal  solution  is 


X  =  [6,  0,  1,  1,  1,  0,  1,  1,  0]  . 


VI.  Negative  Coefficients: 

When  negative  coefficients  appear  in  the  last  M  +  L  equations, 
the  algorithm  is  changed  in  a  quite  obvious  way.  We  can  assume  without 
loss  of  generality  that  the  coefficients  are  +1  or  -1  and  the  last 
L  right  hand  side  components  are  +1.  Theorems  1,  2,  and  3  still  hold. 
For  each  inessential  set  the  basic  column  must  obviously  have  a  +1  co¬ 
efficient  in  the  last  L  components  for  any  basic  feasible  solution. 

In  the  pricing  process  if  the  column  A^  to  be  priced  has  a  negative 
coefficient  in  the  last  L  components,  the  appropriate  p.  is  subtracted 
rather  than  added  to  IIA^  .  Also,  if  an  inessential  variable  with  a 
"negative  one"  is  to  be  introduced,  the  set  must  first  be  made  essential. 
Other  than  these  modifications,  the  algorithm  proceeds  exactly  as  before. 
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